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An Analysis of the Shapes of Interstellar Extinction Curves. V. The 
IR-Through-UV Curve Morphology 

E.L. Fitzpatrick 1 , D. Massa 2 
ABSTRACT 

We study the IR-through-UV wavelength dependence of 328 Galactic interstellar extinction 
curves affecting normal, near-main sequence B and late O stars. We derive the curves using a 
new technique which employs stellar atmosphere models in lieu of unreddened "standard" stars. 
Under ideal conditions, this technique is capable of virtually eliminating spectral mismatch errors 
in the curves. In general, it lends itself to a quantitative assessment of the errors and enables a 
rigorous testing of the significance of relationships between various curve parameters, regardless 
of whether their uncertainties are correlated. Analysis of the curves gives the following results: 

1. In accord with our previous findings, the central position of the 2175 A extinction bump is 
mildly variable, its width is highly variable, and the two variations are unrelated. 

2. Strong correlations are found among some extinction properties within the UV region, and 
within the IR region. 

3. With the exception of a few curves with extreme (i.e., large) values of R(V), the UV and 
IR portions of Galactic extinction curves are not correlated with each other. 

4. The large sightline-to-sightline variation seen in our sample implies that any average Galactic 
extinction curve will always reflect the biases of its parent sample. 

5. The use of an average curve to deredden a spectral energy distribution (SED) will result 
in significant errors, and a realistic error budget for the dereddened SED must include the 
observed variance of Galactic curves. 

While the observed large sightline-to-sightline variations, and the lack of correlation among the 
various features of the curves, make it difficult to meaningfully characterize average extinction 
properties, they demonstrate that extinction curves respond sensitively to local conditions. Thus, 
each curve contains potentially unique information about the grains along its sightline. 

Subject headings: ISM:dust, extinction — methods:data analysis 



1. INTRODUCTION 

In the previous paper in this series (Fitzpatrick & Massa 2005a, hereafter Paper IV), we introduced a 
technique, "extinction- without-standards," to determine the shapes of UV-through-IR interstellar extinction 
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curves by modeling the observed spectral energy distributions (SEDs) of reddened early-type stars. The 
method involves a ^-minimization procedure to determine simultaneously the basic properties of a reddened 
star (namely, T e g, logg, [m/H], and Vturb) o-nd the shape of its extinction curve, utilizing grids of stellar 
atmosphere models to represent intrinsic SEDs and an analytical form of the extinction curve, whose shape 
is determined by a set of adjustable parameters. 

In general, the benefits of extinction-without-standards are increased accuracy and precision (in most 
applications) over results generated using the standard Pair Method of extinction curve determination and, 
very importantly, a reliable estimate of the uncertainties in the resultant extinction curves. Specifically, the 
advantages of the new method include: (1) the elimination of the requirement for observations of unreddened 
spectral standard stars, (2) the near-elimination of "spectral mismatch" as a source of extinction curve error, 

(3) the ability to produce accurate curves for much more lightly-reddened sightlines than heretofore possible, 

(4) the ability to derive accurate ultraviolet curves for later spectral types (i.e., to late-B classes) than 
previously possible, and (5) the ability to provide quantified estimates of the degree of correlation between 
various morphological features of the curves. 

The chief limitation of the method is that the intrinsic SEDs of the reddened stars must be well- 
represented by model atmosphere calculations. This currently eliminates from consideration such objects 
as high luminosity stars, Be stars, Wolf-Rayet stars, and any spectrally peculiar star. This restriction also 
affects the Pair Method, since there are only a small number of unreddened members of these classes which 
can serve as "standard stars," and it is difficult to be certain that the intrinsic SEDs of the standard stars 
really represent those of the reddened stars. In addition, the extinction-without-standards technique requires 
well-calibrated (and absolutely-calibrated) SED observations. 

In this paper, we apply our extinction-without-standards technique to a sample of 328 Galactic stars for 
which multi-wavelength SEDs are available. For these stars, we derive normalized UV-through-IR extinction 
curves, sets of parameters which describe the shapes of the curves, and sets of parameters which characterize 
the stars themselves. The scope of our discussion focuses on two objectives: (1) the presentation of the 
broad-ranging results and a description of the methodology employed, and (2) a thorough examination of 
general extinction curve morphology. Among the issues addressed is the correlation between the IR and 
UV properties of extinction (see, e.g., Cardelli, Clayton, & Mathis 1989). This paper is not intended as 
a review of Galactic extinction and we explicitly restrict our attention only to issues touched upon by our 
new analysis. The results, however, do constitute a broad view of Galactic extinction and lend themselves 
to numerous other investigations, including the general properties of extinction curves, regional trends in 
extinction properties, the correlation between extinction and other interstellar properties, determination of 
intrinsic SEDs for objects in clusters containing survey stars, the study of small scale spatial variations in 
dust grain populations from stars in cluster extinction curves, etc. We plan to pursue some of these in future 
papers. Preliminary results from an early version of this study were reported by Fitzpatrick (2004; hereafter 
F04), which also provided a review of then-recent progress in interstellar extinction studies. 

The sample of stars chosen for this survey and the data used are described in $2] This is followed in f}3] 
by a brief description of the extinction- without-curves technique, including a discussion of the error analysis 
— which is critical to the analyses in the latter parts of the paper. The results of the survey, essentially 
an atlas of Galactic extinction curves, is presented in <QJ with the full sets of tables and the full set of 
figures from this section available in the electronic edition of the Journal. In Sj5]we briefly discuss the stellar 
properties derived from our analysis and then in $6] present a detailed description of Galactic extinction 
curve morphology, from the IR to the UV spectral regions. Finally, in ij7]we provide a brief summary of the 
chief conclusions of this study. 
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2. THE SURVEY STARS AND THEIR DATA 



In principle, the extinction-without-standards technique can be applied to, or expanded to include, any 
type of SED data. In practice, however, we have developed the technique with specific datasets in mind, 
namely, low-resolution UV spectrophotometry from the International Ultraviolet Explorer (IUE) satellite and 
ground-based optical and near-IR photometry. In this particular study, we utilize IUE spectrophotometry, 
UBV photometry, and JHK photometry from the Two-Micron All Sky Survey (2MASS). Although some of 
the survey stars have additional data available, e.g., Stromgren uvbyfi photometry, we elected — for the 
sake of uniformity — to include only the UBV data, in the analysis. As a result, the stellar parameters we 
derive in this paper will be less accurate than those determined in Paper IV. Nevertheless, the errors in these 
parameters are still well-determined. 

The most restrictive of the datasets is that of the IUE, and so we began our search for survey stars with 
the IUE database. Using the search engine provided by the Multimission Archive at STScI (MAST), we 
examined all low- resolution spectra for stars with IUE Object Class numbers of 12 (main sequence O), 20 
(B0-B2 V-IV), 21 (B3-B5 V-IV), 22 (B6-B9.5 V-IV), 23 (B0-B2 III-I), 24 (B3-B5 III-I), and 25 (B6-B9.5 III- 
I). Because our goal is to obtain a uniform dataset of reddened stars whose SEDs can be modeled accurately, 
we eliminated the following types of objects from the available field of candidates: (1) stars without good- 
quality spectra from both the short-wavelength (SWP) and long-wavelength (LWR or LWP) IUE spectral 
regions, (2) clearly-unreddened stars, (3) known Be stars, (4) luminosity class I stars, (5) O stars more 
luminous than class V, (6) O stars earlier than spectral type 05, and (7) stars with peculiar-looking UV 
spectra (as based on our own assessment). 

The list of potential candidates from the IUE database was then examined for the availability of UBV 
data, using the General Catalog of Photometric Data (GCPD) maintained at the University of Geneva 
(see Mermilliod, Mermilliod, & Hauck 1997)0. In general, stars without V and B — V measurements were 
eliminated from consideration. However, broadband Geneva photometry was available for five stars without 
UBV data (BD+56 526, HD62542, HD108927, HD110336, and HD143054), and so we included these stars 
and utilized the Geneva U — B and V — B indices. 

The final trimming of the survey sample was not performed until the SEDs had been modeled. At this 
point we imposed the requirement that all survey stars must have values of E(B — V) > 0.20 mag. This 
limit is somewhat arbitrary and, as shown by Paper IV, useful extinction curves can be derived via the 
extinction-without-standards technique for E(B — V) values considerably lower than 0.20 mag. However, 
the uncertainties do rise at low E(B — V) and for this survey we wanted a sample of stars for which the 
uncertainties in the final parameters were uniformly small. We plan in the future to examine the extinction 
properties along lightly reddened sightlincs. 

Near-IR JHK photometry (and their associated uncertainties) were retrieved for the survey stars from 
the 2MASS database at the NASA/IPAC Infrared Science Archive (IRSA)El Only those 2MASS magnitudes 
for which the uncertainties are less than ±0.1 mag are used here. For seven stars whose 2MASS measurements 
were either non-existent or of low quality (HD23180, HD23512, HD37022, HD144470, HD147165, HD147933, 
and HD149757), Johnson JHK magnitudes were available and were retrieved via the GCPD catalog. Note 
that the availability of JHK data was not a requirement for inclusion in this survey although, as will be 



1 The GCPD catalog was accessed via the website at http://obswww.unige.ch/gcpd/gcpd.html 
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quantified below, most of our sample have such data. 

The last data collection activity was to retrieve ancillary information for all stars, such as coordinates, 
alternate names, and spectral types, using the SIMBAD database. 

The only data processing required for this program involved the IUE spectrophotometry. The IUE data 
contained in the MAST archive were processed using the NEWSIPS software (Nichols & Linsky 1996). As 
discussed in detail by Massa and Fitzpatrick (2000; hereafter MF00), these data contain significant thermal 
and temporal dependencies and suffer from an incorrect absolute calibration. We corrected the data for 
their systematic errors and placed them onto the HST/FOS flux scale of Bohlin (1996) using the corrections 
and algorithms described and provided by MF00. This step is absolutely essential for our program since 
our "comparison stars" for deriving extinction curves are stellar atmosphere models and systematic errors in 
the absolute calibration of the data do not cancel out as they would in the case of the Pair Method. (Note 
that the thermally- and temporally-dependent errors in the NEWSIPS data would not generally cancel out 
in the Pair Method — see MF00.) When multiple spectra were available in one of IDE's wavelength ranges 
(SWP or LWR and LWP), they were combined using the NEWSIPS error arrays as weights. Small aperture 
data were scaled to the large aperture data and both trailed and point source data were included. Short 
and long wavelength data were joined at 1978 A to form a complete spectrum covering the wavelength range 
1150 < A < 3000 A. Data Ion gward of 3000 A were ignored because they are typically of low quality and 
subject to residual systematic effects. 

After all the limits and restrictions were imposed, we arrived at our final sample of 328 stars for which 
UV spectrophotometry covering 1150-3000 A and UB V photometry are available. Of these, 298 stars have 
at least some near-IR photometry, while 287 have a complete set of J-, H-, and K- band measurements 
(with 280 of these from the 2MASS program) . Table Q] lists all the survey stars, along with some general 
descriptive information. (The complete version of Table [T] appears in the electronic version of the paper.) 
The stars in Table 1 are ordered by Right Ascension. For the star names, we adopted the most common 
form among all the possibilities listed in the SIMBAD database (i.e., "HDnnn" was the most preferred, 
followed by "BDnnn" , etc.) There are 185 survey stars which are members of open clusters or associations. 
The identity of the cluster or association is cither contained in the star name itself (e.g., NGC 457 Pesch 
34) or is given in parentheses after the star's name. The spectral types in Table [T] were selected from those 
given in the SIMBAD database, and the source of the adopted types is shown in the "Reference" column 
of the table. When multiple types were available for a particular star, we selected one based on our own 
preferred ranking of the sources. For the B stars, the quality of the spectral types varies widely, and the 
types themselves are given only as a general reference — they do not play any role in our analysis. A scan 
of the types seems to indicate that in some instances we have violated our selection criteria, e.g., several 
"Ib"s and "e"s can be found. However, in our estimation, these are unreliable types and do not reflect the 
spectral information we have examined. For example, a number of B stars in clusters have been erroneously 
classified as emission- line stars based on contamination of their spectra by nebular emission lines. The O 
stars in the sample are, on the other hand, uniformly well-classified and the types are used in the analysis 
(see SI). 

An overview of the survey sample can be gained from Figures [1] and [51 The former shows the location 
of the stars on the sky (in Galactic coordinates) and the latter summarizes the breadth of the stellar and 
interstellar properties of the sample. Figure [T] shows that our sample is clearly biased to sightlines passing 
through the galactic plane, as would be expected given by the lower limit of E(B — V) = 0.20 mag which we 
imposed. The locations of open clusters or associations for which five of more members are in the survey 
sample are indicated on the figure by the large circles. Note that the size of the circles is not intended to 
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represent the physical extent of the clusters or associations. The data shown in histogram form in Figure 
[5] are final results from the analysis, but are useful here to characterize the sample. Most of our stars are 
mid-to-early B stars (T cS = 15000-30000 K) with a median reddening of E(B-V) = 0.45 mag. The median 
value of R(V) (= Ay/E(B — V)) for the sample is 3.05, essentially identical to the Galactic mean value for 
the diffuse ISM, and our sample is dominated by sightlincs though the diffuse ISM. However — due to the 
relatively small sample size and the biases present in the IUE satellite's choice of targets over the years — 
our survey does not necessarily constitute a representative sample of the various types of regions present in 
the ISM. Care must be taken in interpreting average properties derived from our results. 

3. THE ANALYSIS 

It was shown by Paper IV that the energy distributions of reddened B-type stars could be modeled 
successfully using theoretical predictions of the intrinsic SEDs of the stars and a parametrized form of 
the UV-through-IR extinction curve to account for the distortions introduced by interstellar extinction. A 
byproduct of the fit is a determination of the wavelength-dependence of the extinction affecting the star. 
This is the essence of the extinction-without-standards technique. Although it was discussed in detail by 
Paper IV, for completeness, we use this section to outline the basics of the technique. In addition, some 
details of the process have changed since Paper IV, as a result of experience gained from the application of 
the process to several hundred reddened stars. These changes are highlighted here. 

3.1. Modeling the SEDs 

The observed SED f\ of a reddened star can be represented as 

h = Fx eliQ-w B - v ^ x - v ^ R ^ , (i) 

where F\ is the intrinsic stellar surface flux, 0r = R/d is the stellar angular radius (where R is the physical 
radius and d is the distance), E(B — V) is the familiar measure of the amount of interstellar reddening, 
k(X - V) = E(\ - V)/E(B - V) is the normalized extinction curve, and R(V) = A X /E(B - V) is the 
ratio of reddening to extinction at V . By adopting stellar atmosphere models to represent F\ and a using 
parametrized form of k{\ — V), we can treat Equation (fT]) as a non-linear least squares problem and solve for 
the set of optimal parameters which generate the best fit to the observed flux. As in Paper IV, we perform the 
least squares minimization using the Interactive Data Language (IDL) procedure MPFIT developed by Craig 
Markwardt The observed SEDs which are fitted in the process consist of the IUE UV spectrophotometric 
fluxes, optical UBV magnitudes, and near-IR JHK magnitudes discussed 

As related in Paper IV, we developed this analysis utilizing R.L. Kurucz's (1991) line-blanketed, hy- 
drostatic, LTE, plane-parallel ATLAS9 models, computed in units of erg cm -2 sec -1 A -1 and the synthetic 
photometry derived from the models by Fitzpatrick & Massa (2005b). The models are functions of four 
parameters: T c ff, logg, [m/H], and v%. All of these parameters can be determined in the fitting process 
although, because of data quality, it is sometimes necessary to constrain one or more to a reasonable value 
and solve for the others. We smooth and bin the IUE fluxes to match the sampling of the ATLAS9 models 
(10 A bins over most of the IUE range; see Fitzpatrick & Massa 1999). 



3 Markwardt IDL Library available at http://astrog.physics.wisc.edu/~craigm/idl/idl.html 
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Because we include some O stars in the current sample, we have expanded our technique to incorporate 
the TLUSTY OSTAR2002 grid of line-blanketed, hydrostatic, NLTE, plane-parallel models by Lanz & 
Hubeny (2003). That grid includes 12 T c g values in the range 27500 - 55000 K, 10 chemical compositions 
from twice-solar to metal-free, and surface gravities ranging from log g = 4.75 down to the modified Eddington 
limit. All models were computed with v tur b = 10 kms -1 . We only consider the solar abundance models in 
this analysis and, thus, the TLUSTY models are considered as functions of two parameters, T B s and log g. 
Synthetic UBV photometry for these models was produced as described above for the ATLAS9 models. 
To keep the O-star and B-star fitting procedures as similar as possible, the TLUSTY models were binned 
to the same wavelength scale as the ATLAS9 models. While analyzing the O stars we found that their 
low-dispersion IUE spectra (excluding the strong wind lines) are very insensitive to temperature and our 
analysis yielded very uncertain results. As a result, we modified the procedure for these stars and adopted 
values of T c g based on their spectral types, rather than solving for T e g. Table [2] lists the temperature scale 
used; it is a compromise between the results of Martins, Schaerer, & Hillier (2005) and our own analysis 
of optical line spectra of O-stars (such as published by Walborn & Fitzpatrick 1990) using the TLUSTY 
models, which are appropriate for O stars without massive winds. The results of this investigation will be 
reported elsewhere. We assume a generous uncertainty in these T e g values (see below) and the extinction 
curve results are actually very insensitive to the adopted temperatures. 

The value of the surface gravity log g is often poorly-determined when using only broadband photometry, 
because it lacks a specific gravity-sensitive index to help constrain logg. For the cluster stars, however, we 
can apply ancillary information — specifically, the cluster distances as listed in Table Q] — to provide strong 
constraints on \ogg. We adopt the same procedure as used by Fitzpatrick & Massa (2005b), in which the 
Padova grid of stellar structure models allow the Newtonian gravity of a star to be inferred through its 
unique relation with the star's surface temperature and radius. (When the distance is specified, the physical 
radius of the star becomes a fit parameter via its influence on the angular radius 6r. Fitzpatrick and Massa 
(2005b) used distances determined by Hipparcos parallaxes.) In our iterative fitting procedure, the current 
values of T e // and Or, coupled with the Padova models, determine the current value of \ogg. Generous 
\-<j uncertainties in the distances are included in the error analysis (see below). For field stars, no such 
constraints on log g are possible and we solve for logg as a free parameter. We do, however, apply a reality 
check to the results and, if the final value seems physically unlikely (i.e., \ogg > ~4.3 or logg < ~3.0), we 
replace it with the mean sample value of logg = 3.9. Uncertainties in this assumed value are incorporated 
in the error analysis. 

A parametrized representation of the extinction curve, covering the whole UV-through-IR spectral range, 
is the heart of the current analysis. As in Paper IV, we construct this curve in two parts, joined together at 
2700 A. An example of this formulation, for one particular set of parameters, is illustrated in Figure [3] In 
the UV, (A < 2700 A; the shaded region in the figure), we use a modified form of the UV parametrization 
scheme of Fitzpatrick & Massa (1990; hereafter Paper III), as shown by the thick solid curve. At longer 
wavelengths, we use a cubic spline interpolation through a set of UV, optical, and IR "anchor points" (the 
U n , O ni and I n in the figure), as shown by the thick dashed curve. 

The modified Paper III extinction curve is defined by 



where x = A 1 , in units of inverse microns (/im 1 ) . There are seven free parameters in the formula which 
correspond to three features in the curve: (1) a linear component underlying the entire UV wavelength range, 




(2) 
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defined by c\ and C2; (2) a Lorentzian-like 2175 A bump, defined by C3, xq, and 7 and expressed as 



(3) 



and (3) a far-UV curvature component (i.e., the departure in the far-UV from the extrapolated bump-plus- 
linear components), defined by C4 and C5. All seven free parameters can be determined by the least squares 
minimization algorithm. 

The modification made here to the Paper III formula is in the far-UV curvature term. In Paper III, the 
value of C5 was fixed at 5.9 /nm" 1 and C4 was the scale factor applied to a pre-defined cubic polynomial. In 
working with the current dataset, we found that the formulation in Equation @ significantly improved the 
fits to many stars, particularly those with weak far-UV curvature, and degraded the fits in almost no cases. 
In the modified form, the curvature is functionally simpler — containing a single quadratic term — although 
we have added another free parameter to the extinction curve representation. Because the primary goal of 
the Paper III formula was (and still is) to provide an analytical expression which reproduces as closely as 
possible observed extinction curves, the cost of an additional free parameter was deemed worthwhile. 

Using the UV fit parameters above, additional quantities can be defined which help describe the UV 
curve properties. Particularly useful ones include 1) A1250 = 04(8.0 — C5) 2 , which is the value of the FUV 
curvature term at 1250 A and provides a measure of the strength of the FUV curvature; 2) Ab ump = n 03/(27), 
which is the area of the 2175 A bump; and 3) Eb ump = C3/7 2 , which is the maximum height of the 2175 A 
bump above the background linear extinction. 

The cubic spline interpolation which produces the optical-through-IR region of our parametrized ex- 
tinction curve is produced using the IDL procedure SPLINE. The nine anchor points shown in Figure [3] are 
specified by five free parameters. The UV points, U\ and U%, are simply the values at 2600 A and 2700 
A resulting from the modified Paper III formula and require no new free parameters. They assure that 
the two separate pieces of the extinction curve will join smoothly, although not formally continuously. The 
three optical points, 0\, O2, and O3 (at 3300 A, 4000 A, and 5530 A, respectively) are each treated as free 
parameters and are adjusted in the fitting procedure to assure the normalization of the final extinction curve. 
The IR points, 1%, I2, I3, I4, and I5 (at 0.0, 0.25, 0.50, 0.75, and 1.0 /um -1 , respectively) are functions of 
two free parameters, km and R(V) as follows: 



This assures that the IR portion of our curve follows the the power-law form usually attributed to IR 
extinction, with a value for its exponent from Martin & Whittet (1990). As noted in Paper IV, the value of 
the power law exponent could potentially be determined from an analysis like ours. However, we have found 
that an IR dataset consisting only of JHK magnitudes, as in our survey, is insufficient to specify three IR 
parameters. 

Note that we adopt the cubic spline formulation for the optical/IR extinction curve simply because we 
do not have an acceptable analytical expression for the curve shape over this range. The spline approach is 
very flexible in that the number of anchor points can be modified depending on the datasets available. In the 
current application, we use only three optical anchor points because we have only three optical data points 
(V, B — V, and U — B). This approach will bear its best fruit when it can be applied to spectrophotometric 
data in the near-IR through near-UV region. Then, a large number of anchor points can be used to precisely 
measure the heretofore poorly-determined shape of extinction in this region, without bias towards a particular 
analytical expression. 



7„ = k(X -V) = k m X 



-1.84 



R(V) 



(4) 
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In summary, the analysis performed here — modeling the SEDs of 328 reddened stars via Equation ((T|) 
— involves determining the best-fit values for as many as 18 free parameters per star, via a non-linear least 
squares analysis. These include up to 4 parameters to define the theoretical stellar atmosphere model (T e a, 
logg, [m/H], Vturb), up to 12 parameters to describe the extinction curve shape (0\, O2, O3, R{V), km, 
c\ through C5, xq, 7), the angular radius 8r, and E(B — V). We weight the UV, optical, and IR datasets 
equally in the fitting procedure. 

3.2. Error Analysis 

One of the main benefits of the extinction-without-standards technique is the error analysis, which 
provides a well-quantified estimate of the uncertainties in the best-fit model parameters and allows possible 
correlations between parameter errors to be explored. This latter benefit is important for assessing the reality 
of apparent correlations between parameters. The uncertainties in the best-fit parameters derived for our 
survey stars were determined by running 100 Monte Carlo simulations for each star. In each simulation, the 
fitting procedure was applied to an input SED consisting of the final best-fit model convolved with a random 
realization of the observational errors expected to affect the actual data. The adopted l-a uncertainties for 
each parameter, which will be presented in 21 were taken as the standard deviations of the values produced 
by the 100 simulations. 

Our observational error model for the IUE data consists of random photometric uncertainties and 
camera zero-point errors as described by Fitzpatrick & Massa (2005b). The assumed observational errors 
in the Johnson B — V and U — B indices were as given in Table 7 of that paper. The V magnitudes were 
assumed to have a l-a uncertainty of 0.015 mag. The uncertainties in the 2MASS JHK data were as 
obtained from the 2MASS archive. Johnson JHK magnitudes were assumed to have uncertainties of ±0.03 
mag. Random realizations of each of these observational errors, which were added to the best-fit model SED 
for each Monte Carlo simulation, were determined using the IDL procedure RANDOMN, which produces a 
normally distributed random variable. 

In cases where assumptions were made about the values of specific fit parameters, we incorporated 
uncertainties in the assumptions in the error analysis. In particular: (1) for the O-type stars, the adopted 
spectral type-dependent temperatures were taken to have l-a uncertainties of ^1000 K; (2) for cluster stars, 
the adopted distances (used to constrain the surface gravities) were assumed to have l-a uncertainties of 
±20%; and 3) for the field stars whose values of \ogg were taken to be the sample mean of 3.9, this mean 
was assumed to have a l-a uncertainty of ±0.2. The values of T e g, distance, and log g used in the the Monte 
Carlo calculations for the relevant stars were varied randomly in the simulations (using RANDOMN) in 
accord with these uncertainties. 

This study certainly does not constitute the first attempt to quantify the uncertainties in interstellar 
extinction curves. Most pair method studies (see, for example, Cardelli, Sembach, & Mathis 1992) have 
incorporated some form of error analysis, often based on the methodologies presented by Massa, Savage, 
& Fitzpatrick (1983) and Massa & Fitzpatrick (1986). However, as long as we restrict our sample to stars 
which are well represented by the model atmospheres we employ, the advantages of the current technique are 
great. Because the stellar parameters (temperature, surface gravity, and abundance) are given by continuous 
mathematical variables (instead of a non- uniformly sampled, discrete sets of standard stars), we are able to 
perform a well-defined Monte Carlo analysis. The results of this analysis explicitly quantify the uncertainties 
in all of the input data and assumptions and, thus, the final error bars affecting the derived curves. Moreover, 
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since many realizations of the individual curves are produced, the full shape of the "error ellipses" (which 
describe correlations between the errors) are determined for each specific set of input parameters. Additional 
discussion of the extinction- without-standards error analysis can be found in Paper IV of this series, along 
with a demonstration of the quantitative accuracy of the results. 

4. AN ATLAS OF GALACTIC EXTINCTION CURVES 

The results of the extinction-without-standards analysis of 328 Galactic stars arc presented in Table 
[21 Tabled) an d Figure dl The 18 free parameters determined by the fitting procedure are divided between 
Tables [3] and d] with the latter containing the 12 parameters which define the shape of the normalized 
interstellar extinction curves k(X — V) = E(X — V)/E(B — V). For both tables, only the first 10 entries 
are shown here. The full versions of the tables can be viewed in the electronic edition of the Journal. The 
uncertainties listed in the tables are the l-cr errors derived from the Monte Carlo analysis described in §3.21 

Figure 0] shows the normalized extinction curves for the survey sample. The figure consists of 33 panels, 
each (except the last) containing 10 extinction curves arbitrarily shifted vertically for clarity. Only the first 
panel is shown here. The full figure is given in the electronic version of this paper. The solid curves in the 
figure show the parametrized UV-through-IR curves whose shapes were determined by the fitting procedure 
described in §2] (the parameters describing the curves are in Tabled]). An estimate of the shape of the average 
Galactic extinction curve (corresponding to R(V) — 3.1; from Fitzpatrick 1999) is shown for reference by 
the dash-dot curves. The l-cr uncertainties of the survey extinction curves are indicated in Figured] by the 
grey shaded regions, which are based on the Monte Carlo error simulations. Their thicknesses indicate the 
standard deviations of the ensemble of simulations at each wavelength. The actual normalized ratios between 
the observed stellar SEDs and the atmosphere models are shown by the symbols. Large filled circles indicate 
JHK data in the IR region (A -1 < 1 jum -1 ) and UBV data in the optical region (1.5 < A -1 < 3.0 /im -1 ). In 
the UV region ((A -1 > 3.3 /irn^ 1 ), the small symbols with l-cr error bars show the ratios between the IUE 
data and the models. 

A close examination of the curves in Figure [4] shows that the parametrized curves are extremely good 
representations of the observed extinction ratios and thus serve as useful proxies for the actual curves 
themselves. This is particularly apparent in the UV, where the spectrophotometric data show the flexibility 
of the parametrization scheme. For those wishing to use these curves in extinction studies, we have prepared 
a tar file containing the parametrized curves for all 328 stars, sampled at 0.087 /irn^ 1 intervals, and their 
accompanying fit parameters. Directions for retrieving the file are given in the Appendix. 

In using or interpreting these curves it is important to recognize that their shapes in the regions between 
the IR and optical and between the UV and optical are interpolations only and not strongly constrained by 
data. Additional observations, particularly fully-calibrated spectrophotometric data, would be very useful 
to constrain the shape of the extinction curves in these regions — and in the optical where only broadband 
UBV measurements are currently employed. It is certainly counter-intuitive that the spectral regions where 
the detailed shapes of the extinction curves are most poorly-determined are ground-accessible, while the UV 
data are so well-measured. 
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5. Properties of the Sample Stars 

Although the main goal of this paper is to explore Galactic extinction, it is nevertheless reasonable 
to consider briefly the stellar properties revealed by our analysis since they directly affect the extinction 
results. As discussed above, our reliance on broadband photometry for the current work results in stellar 
parameters that are not as accurate as those presented in Paper IV, due to the lack of a good surface gravity 
discriminator. Nevertheless, more than 50% of the sample stars reside in open clusters and associations and, 
for these stars, the accuracy will be increased, due the use of ancillary information. 

Of the four stellar properties determined in the analysis (T c ff, log<?, [m/H], Vturb), the most significant to 
our extinction program is T e e because it has the most impact on the shapes of the model atmosphere SEDs 
used to derive the extinction curves. Figure [5] shows a plot of the derived T e g values against spectral type for 
the spectral class BO and later stars (filled and open circles). For comparison, we show several spectral type 
vs. T c g calibrations from the literature (solid and dashed lines) and data from our photometric calibration 
study (Fitzpatrick & Massa 2005b, open squares) in which we modeled the SEDs of 45 unreddened B-type 
stars. While the survey star data are generally consistent with the comparison data in the figure, considerable 
scatter is present at some spectral classes — most particularly at types B2 and B3 — and there is a general 
departure between our results and the T e g calibrations in the neighborhood of types Bl and B2, in the sense 
that our results indicate hotter temperatures than the calibrations. 

We examined the IUE spectra of a number of the survey stars — those whose temperatures are most 
discrepant with their claimed spectral types — and compared them with spectral classification standards 
to see if somehow our fitting procedure was arriving at grossly incorrect temperatures. An example of such 
a comparison is shown in Figure [6l We plot a portion of the IUE spectra of the survey stars HD228969 
and BD+45 973, the two hottest stars in the "B2" and "B3" spectral bins, respectively, along with several 
spectral standards with expected temperatures in the neighborhood of 30000 K. The close match between 
the survey stars and the hot standards is evident and it is clear that the cool spectral types found in the 
literature for these stars are unreliable. This is the general conclusion from all such comparisons we have 
performed. The temperatures found from the fitting procedure are consistent with those expected from a 
close examination of the UV spectral features of the stars. We conclude that the outliers in Figure [5] result 
from poor optical spectral types. This is not surprising, since the available types are from a large number of 
sources and based on a wide variety of observational material of very non-uniform quality. 

The general discrepancy between our results and the calibrations in the B1-B2 region is a different 
matter. We derive considerably hotter effective temperatures for the Bl (25000 - 26000 K) and B2 (23000 
- 24000 K) stars than expected from previous calibrations. However, inspection of the UV features in the 
ATLAS9 models make it difficult to believe that typical Bl and B2 stars are as cool as the spectral type - 
T e ff calibrations suggest. We must bear in mind, however, that our sample is strongly biased toward cluster 
stars, which may be considerably younger and more compact than the "field" B stars used in the calibrations. 
Furthermore, the current B star calibrations are all over 20 years old, and it is quite possible that they are 
in need of revision. 

Another way to look at the T e s values is shown in Figure where we plot T e g as a function of the 
Stromgren reddening-free index [c] = c\ — 0.20(& — y), which is a measure of the strength of the Balmer 
jump. Stromgren photometry is not used in this program but is available for 162 of our stars. The symbols 
in Figure [7] are the same as in Figure with the addition of the open circles which denote O stars. The 
figure demonstrates the essentially exact overlap between the current results and those for the unreddened, 
mid-to-late B stars from Fitzpatrick & Massa (2005b) as well as the smooth transition from the early B stars 
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into the O stars. There is no indication of any systematic effects present in the results for the early B stars. 
On the contrary, the spectral type vs. T c g relations, transformed into the T Q g vs. [c] plane as described in the 
figure legend, show a number of abrupt and physically unrealistic changes in slope suggestive of inadequacies 
in the calibrations. 

We conclude that our derived effective temperatures are reasonable. In most cases where a temperature 
strongly disagrees with a published MK type, it agrees quite well with the UV type determined by Valencic, 
Clayton, & Gordon (2004), indicating that the MK type is of poor quality or else influenced by something 
else (e.g., the presence of a cooler companion which is invisible in the UV). We also suspect that the spectral 
type - T c g calibration for the Bl and B2 stars may need to be revised. Finally, we are gratified by the overall 
consistency between the current results and those of Paper IV and Fitzpatrick & Massa (2005b) , where the 
stellar parameters were more strongly constrained, and with the smooth relation between T c s and [c], which 
suggests that no strong systematic effects are present. 

6. Properties of Galactic Interstellar Extinction 

6.1. General Properties 

Although the individual extinction curves for all of the survey stars are displayed in the 33 panels of 
Figure [H it is nonetheless difficult — from that figure — to visualize the range of extinction properties 
present in the sample. To provide such a view, we plot the analytical fits for the full set of 328 survey 
curves in Figure [8] (These curves can be reproduced from the parameters listed in Table |4l The curves have 
been plotted using small dots in those spectral ranges where they are interpolated or extrapolated. The 
solid portions correspond to the regions constrained by near-IR JHK photometry, optical UBV photometry, 
and UV IUE spectrophotometry. The top panel, in which the curves are plotted in their native form, 
E(X — V)/E(B — V), shows the wide range of variation observed in Galactic extinction curves, although a 
clear core of much more restricted variation is evident in the distribution. The convergence of the curves in 
the range 1.8 < A -1 < 2.5 /Ltm -1 is a product of the normalization and obscures our view of variations in the 
optical region. The bottom panel of Figure[5]presents the same curves, but normalized to the total extinction 
at 1 fim . The exact overlap of the curves at A" 1 < 1 /im _1 arises because we have adopted a power law 
form with a fixed exponent of —1.84 for all curves in this spectral region (see Eq. [4]). It is unlikely that 
the exponent is actually so constant — Larson & Whittet (2005) have found evidence for a more negative 
value in a sample of high-latitude clouds — however, as noted earlier in 331 our near-IR dataset consisting 
of only JHK measurements is insufficient for independently evaluating the exponent. This IR normalization 
reveals that the extinction in the optical spectral region can range from very-nearly grey to very strongly 
wavelength-dependent. 

One of the goals of many Galactic extinction studies is to derive an estimate of the typical or average 
wavelength-dependence of interstellar extinction. Such a mean curve is often used as the standard of normalcy 
against which particular sightlines are judged, or for comparison with results for external galaxies, or for 
"dereddening" the SEDs of objects for which there is no specific extinction knowledge. While constructing 
an average curve is a straightforward process, the degree to which the result represents "normal" or even 
"typical" extinction is problematic. We will take up the issue of an average Galactic curve later in iJTl Here 
we present an average curve for our sample, to help characterize the general extinction properties of our 
survey sightlines. 

The bottom panel in Figure [2] shows a clear peak in the distribution of R(V) values of our sample. A 
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Gaussian fit to the region of this peak, shown as the smooth curve in the figure, has a centroid at R(V) = 2.99 
and a width given by a — 0.27. The peak is within the range of values considered as average for the diffuse 
ISM (see, e.g., Savage & Mathis 1979). Thus, as noted in our sample is dominated by sightlines whose 
R(V) values are consistent with the diffuse phase of the ISM (although composite sightlines are undoubtedly 
present). We have constructed an average curve to represent the properties of these diffuse sightlines by 
taking the simple mean extinction value at each wavelength using all the curves with 2.4 < R(V) < 3.6 (i.e., 
the 2-tr range of the Gaussian fit in Figure [2 243 sightlines in total). This mean extinction curve is shown in 
Figure [9] by the thick solid curve. The dark grey shaded region shows the variance of the 243-curve sample. 
The set of 12 parameters describing this curve are listed in Table O Since so many of our sightlines are 
included in the mean, removing the restriction on R(V) has little effect. If we had included all 298 sightlines 
for which R(V) has been derived, then the mean curve would differ from that shown only by being several 
tenths lower in the UV region. The variance of the full sample is larger, and this is shown by the lightly 
shaded region. For comparison, several other estimates of average Galactic curves are shown in Figure [HI 
The curves from Cardelli et al. (1989), Fitzpatrick (1999), Seaton (1979), and Valencic et al. (2004) are 
all intended to represent the diffuse ISM mean. The results from Savage et al. (1985) are mean values for 
the 800 sightlines in that study with E(B — V) > 0.20 mag (matching our survey cutoff). No restriction on 
R(V) was imposed. The error bars for the Savage et al. (1985) data are sample variances; they are generally 
similar to the variances of our full sample (lightly shaded region). The much larger value for the 1500 A 
point is likely due to spectral mismatch in the potentially strong C IV stellar wind lines. 

The differences among the various mean curves in Figure [9] are instructive. The great intrinsic variety 
of Galactic extinction curves as seen, for example, in Figure |8] shows that any mean curve is subject to the 
biases in the sample from which it was produced. It is probably impossible to construct a sample of sightlines 
whose properties could be claimed to provide a fair representation of all the types of conditions found in the 
ISM. Thus, there is likely no unique or best estimate of mean Galactic extinction. Any of the mean curves in 
Figure [9] would serve as reasonable representations of Galactic extinction. In any situation where an average 
curve is adopted, however, it is important to recognize the intrinsic variance of the underlying sample and 
incorporate the uncertainties of the average curve in any error analysis. Because it is derived from such a 
large sample and is largely free of contributions from spectral mismatch, the sample variance for our diffuse 
curves (shown in Figure [9] by dark shaded region) would provide a reasonable estimate of the uncertainty 
in any version of a mean Galactic diffuse ISM curve. We have included our diffuse ISM mean curve from 
Figure [9] — and its accompanying uncertainty — in the tar file discussed in §4] (see the Appendix). 

We can also use our large sample to investigate the "smoothness" of UV extinction. In Figure [TUl we 
plot the simple mean of the actual extinction ratios for our survey sample in the spectral region covered 
by the IUE data (upper curve, open circles). Overplotted is a parameterized fit to these data, using the 
extinction formulation given in Sj3] (solid curve). This figure illustrates two points; namely, 1) the lack of 
small-scale structure in UV extinction and 2) the degree to which the modified Paper III UV parametrization 
scheme reproduces the shape of UV extinction. A detailed discussion of the former point, and an indication 
of the kinds of features that might be expected in the UV, can be found in Clayton et al. (2003). In 
short, polycyclic aromatic hydrocarbon (PAH) molecules, which have been suggested as the source of mid- 
IR emission features, might produce noticeable absorption in the UV, and possibly even contribute to the 
2175 A bump absorption. Earlier studies have always failed to find structure in UV extinction curves and 
Clayton et al. were able to place very stringent 3-<r upper limits of ~ 0.02Ay on any possible 20 A-wide 
features in the extinction curves towards two heavily reddened B stars. The data at the bottom of Figure [TOl 
show the differences between the mean extinction curve and the best-fit model. A small number of points 
do rise above (or below) the general noise level of the residuals, but these points — which are labeled in the 
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figure — are due to interstellar gas absorption lines (which are sometimes strong in the observations and 
always non-existent in the model atmosphere SEDs), mismatch between the C IV A1550 stellar wind lines in 
the O stars and the static model SEDs, or known inadequacies in the ATLAS9 opacity distribution functions 
(labeled "b" in the figure; see Paper IV). Excluding these points, the standard deviation of the mean curve 
around the best-fit is 0.06E(B — V) mag, corresponding to ~ 0.02Ay mag, for the 10 A- wide spectral bins. 
This is not quite so restrictive as the result of Clayton et al., but once again affirms the smoothness of UV 
extinction. 



6.2. Spatial Trends 

Studying spatial trends in extinction can serve two purposes. First, identifying strong regional variations 
can allow a better estimate of the shape of the extinction curve affecting a particular sightline permitting, for 
example, a more accurate determination of the intrinsic SED shape of an exotic object. This has been the 
motivation for a number of studies, such as by Massa & Savage (1981) and Torres (1987) who used extinction 
curves derived for B stars in open clusters (NGC 2244 and NGC 6530, respectively) to determine the SEDs 
of cluster O stars. The second purpose is to gain insight into the nature of dust grains and the processes 
which modify them by observing how extinction curves respond to various environmental properties, such 
as density or radiation field. 

In Figure [TT] we show a sweeping view of the regional trends in our data by plotting several extinction 
curves properties against Galactic longitude for each of our survey sightlines. Sightlines towards stars in the 
clusters or associations highlighted in Figure Q] are indicated by the larger symbols in Figure [IT] (the key is 
in figure caption) and the rest of the sample by the small filled circles. The dashed lines in the figure panels 
show the parameter values that correspond to the diffuse mean curve in Figure [9] 

A number of regions stand out in Figure [TTJ For example, NGC 1977 near I = 210° (open triangles; 
includes the Orion Trapezium region) is well-known for its high R(V) sightlines. Figure [TT1 shows that this 
region is also notable for its flat UV extinction curves (i.e., small C2) and weak bumps (i.e., small Ab ump ). 
Interestingly, the far-UV curvature for this region appears typical. An early discussion of UV extinction 
curves towards NGC 1977 can be found in Panek (1983). The relationship between R(V) and other curve 
properties will be the subject of £16.41 below. 

The Carina direction (large filled circles), particularly towards Tr 14 and Tr 16, also shows elevated 
R(V) values. This is another weak-bumped direction which also shows lower than average FUV curvature. 
Optical and IR extinction studies have been performed for Carina sightlines (e.g., Tapia et al. 1988 and 
references within), but we are not aware of correspondingly detailed UV extinction studies. 

Large R(V) values are also seen along a number of sightlines in the general direction of the Galactic 
center. This includes sightlines to the cluster NGC 6530 (open squares) and towards the p Oph dark cloud 
(small filled circles near I ~ 253°. These sightlines also feature low UV extinction and, in the case of NGC 
6530, weaker than average 2175 A bumps. Large R(V) are well known in the Ophiuchus region (Chini & 
Krugel 1983) and the UV extinction has been examined by Wu, Gilra, & van Duinen (1980). UV extinction 
towards NGC 6530 was studied by Torres (1987). 

Finally we note a broad region from I ~ 50° to I ~ 150° where the R(V) values are systematically 
slightly below the mean value. No trends are obvious in the UV parameters. These sightlines are in the 
direction of the Perseus spiral arm (e.g., Georgelin & Georgelin 1976) and sample dust in the interarm region 
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and, possibly, the Perseus Arm itself — for those stars more distant than ^2 kpc, such as in h & x Per (x's 
in Figure [TTj) . Extinction towards individual regions in this zone have been studied (e.g., Tr 37 by Clayton 
& Fitzpatrick 1987, Cep OB3 by Massa & Savage 1984, h & \ Per by Morgan, McLachlan, & Nandy 1982, 
and Cr 457 by Rosenzweig & Morrison 1986), but we are not aware of any comprehensive investigation of 
the general region. 

Morgan et al. (1982) found a dependence of UV extinction on Galactic latitude, 6, for sightlines to stars 
in h & x Per, in the sense that the extinction at 1550 A increases with increasingly negative values of b. Our 
data suggest a similar effect for the UV linear slope C2 (which is closely related to the extinction level at 1550 
A), but not for any other extinction parameter, including R(V). We have examined our dataset for general 
trends with Galactic latitude, or with distance above and below the plane, and have found none. This is 
not surprising, however, since our sample is dominated by low-latitude sightlines and we have little leverage 
for a latitude-dependence study. Our lower cutoff of 0.20 mag in E(B — V) eliminated most high-latitude 
sightlines from consideration. Previous studies have shown that it is difficult to uncover latitude dependences 
in extinction (e.g., Kiszkurno-Koziej & Lequeux 1987). Local trends might be uncovered by examining small 
zones in Galactic longitude (such as in the study by Morgan et al.), although such detailed investigations 
are beyond the scope of this paper. Likewise, studies of extinction variations over small spatial scales, such 
as among sightlines to cluster members are beyond our scope, but might well provide important information 
linking grain populations with other ISM diagnostics. 

The quantity E(B — V)/d, which can be computed from the data in Tables [T] and [3j provides a crude 
measure of one important physical property of the ISM, namely, dust density. While the shortcomings of 
this measure as a direct proxy for density are clear — for example, a high density dust cloud along a long, 
otherwise vacant sightline will yield a misleadingly low value of E(B — V)/d — it is nevertheless useful as a 
first-look and as a guide for future studies. Figure [T2l shows plots of four extinction curve parameters against 
E(B — V)/d. The symbols are the same as for Figure [TT] The three UV parameters all show evidence for a 
weak trend with density, in the sense of flatter slopes, broader bumps, and increasing FUV curvature with 
increasing density. Hints of these three effects were seen in the first two papers in this series (Fitzpatrick & 
Massa 1986, 1988; hereafter Papers I and II). We see no evidence for a trend with bump strength Af, ump (not 
shown in Figure fT2|) and the dependence of R(V) on E(B — V)/d is complex and difficult to characterize, 
although the highest density sightlines all have larger-than-average values of R(V). While not conclusive, 
the results in Figure [12] certainly suggest that comparisons of our survey data with detailed measures of 
ISM physical conditions could yield interesting results. We note that Rachford et al. (2002) have found 
positive trends of bump width and far-UV curvature with the fraction of hydrogen atoms in the form H2. 
This is consistent with the results in Figure [12] in the sense that one would expect higher H2 fractions in 
denser, and therefore better-shielded, regions. Other studies of the density-dependence of extinction have 
been performed by Massa (1987) and Clayton, Gordon, & Wolff (2000). 



6.3. Relationships Among the Fit Parameters 

With twelve parameters to describe the UV-through-IR extinction curves of each of the 328 survey 
stars, there are many possible correlations and relationships to investigate. We have looked at all of these 
possibilities and, as an interested reader can verify from the data in Table [4] virtually all the parameters 
are remarkably UNcorrelated with each other! In this section, we consider only the two most striking 
relationships between extinction parameters: c\ vs. C2 (the UV linear intercept and slope) and R(V) vs. 
km (the ratio of selective-to-total extinction and the IR scale factor). We also examine the most important 
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non-relationship between parameters: xo vs. 7 (the centroid and width of the 2175 A bump). Below, in 
i )6.4[ we will consider the possible relationship between IR curve features and UV features. 

6.3.1. c\ vs. C2 

In Papers I and II we showed that pair method extinction curves in the region of the 2175 A extinction 
bump could be modeled very precisely using a Lorentzian-like "Drude profile" (see Eq. [3] ) combined with a 
linear background extinction (defined by the parameters c\ and ci in Eq. [2]). Further, it was shown that the 
linear parameters for the 45 sightlines in the study appeared to be very well correlated, and could likely be 
replaced with a single parameter without loss of accuracy. In the pilot study for this paper, F04 showed that 
this correlation was maintained in a larger sample of 96 curves derived using the extinction-without-standards 
technique. 

Figure [T3] shows a plot of the linear parameters C\ vs. c-i for our survey sample of 328 extinction-without- 
standards curves. The dotted error bars show the orientation of the 1-a error ellipses of the measurements. 
These were determined by the distribution of results from the Monte Carlo simulations. The obvious cor- 
relation between the errors in C\ and C2 is not unanticipated (see Figure 4 in Paper II), but the ability to 
explicitly determine such errors is a major advantage of the extinction-without-standards approach and is 
critical for evaluating the significance of apparent correlations. 

The solid line in Figure [13] corresponds to the linear relation 

ci = 2.09 - 2.84 c 2 , (5) 

which is a weighted fit that minimizes the scatter in the direction perpendicular to the fit. (Because there is 
uncertainty in both c\ and C2, a normal least squares fit is not appropriate.) The fact that this relationship 
is nearly parallel to the long axis of the correlated error bars explains why the relationship between c\ and 
C2 remains so clear, even in the presence of observational error. 

To determine whether the observed scatter about the mean ci vs C2 relationship is caused only by 
observational error or is at least partially the result of "cosmic" scatter, we must examine the residuals to 
the best-fit relationship. The distribution of residuals perpendicular to the best-fit line is complex, consisting 
of a Gaussian core of values with a ~ 0.07 and a more extended distribution of outliers reaching out to values 
of about ±0.3. About 86% of the points fall within the 2-a range of the Gaussian core. The RMS value of 
the expected observational errors perpendicular to the best-fit relation is 0.057. The Gaussian core of the 
observed residuals is thus only slightly broader than that expected from observational scatter alone. We 
conclude that c\ and C2 are indeed intrinsically well-correlated quantities, with a cosmic scatter comparable 
to our measurements errors. However a significant fraction of the sample (~10%) show evidence for a wider 
deviation from the mean relationship. In most instances, Equation ([2]) could be simplified without loss of 
accuracy by replacing the two linear parameters ci and C2 with a single parameter. 

6.3.2. R(V) vs. k IR 

F04 showed that the two parameters describing the IR portion of the extinction-without-standards 
curves, km and R(V) are apparently well-correlated and that the IR curve might be defined by a single 
parameter. Figure fl4l shows a plot of fc/_R vs. R(V) for the survey sample, along with the 1-a error bars. 
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The solid curve shows the best-fit weighted linear relationship 



km 



0.83 + 0.63 R(V) , 



(6) 



which minimizes the scatter in the direction perpendicular to the relation. As in the discussion above, we 
see that the errors are strongly correlated and the near coincidence of the long axis of the error ellipses with 
the direction of the best-fit relationship would preserve the appearance of a correlation even in the face of 
significant observational error. The residuals in the direction perpendicular to the best-fit line are distributed 
in a Gaussian form, with a = 0.11. The RMS value of the observational errors in this same direction is 
±0.12. Thus the observed distribution of points in Figure[l4]is consistent with perfectly correlated quantities 
and the expected observational error. 

The correlation between R{V) and km indicates that the shape of near-IR extinction at A > ~1 /xm, 
over a wide range of R{V) values, can be characterized by a single parameter. I.e., the two parameters in 
Equation Q are redundant and could be replaced — without loss of accuracy — by a single parameter, e.g., 
R(V), based on Equation ©: 



This is consistent with the results of Martin & Whittet (1990; see their Table 2), who utilized IR data out 
to the M band near 5 /im. Our study utilizes only shorter wavelength near-IR bands, but includes many 
more sightlines and spans a wider range in R(V) values than could be studied by Martin and Whittet. Our 
parameter km is essentially the same as their e parameter, and our results demonstrate the dependence of 
e on R(V). 

The significance of the tight correlation between R(V) and km is somewhat difficult to assess. We 
must keep in mind that, at one level, this relation simply states that the three IR data points given by 
2MASS JHK photometry can be summarized by two parameters at the level of the errors in the 2MASS 
photometry. On the other hand, it also bears on the issue of the underlying shape of infrared extinction and 
its so-called universality. Our data, which are dominated by diffuse ISM sightlines, are consistent with the 
notion of a universal shape for IR extinction, as given by Equation ([7]). However, given our restricted IR 
wavelength coverage and the typical uncertainties in the 2MASS data, our results are relatively insensitive to 
departures from universality. For example, the actual IR power law exponent among our sample stars could 
vary significantly around the mean value of —1.84 used here, and we would still find a very strong correlation 
between km and R(V). The issue of the universality of IR extinction is best left to indepth studies which 
utilize more focussed approaches and more appropriate datasets, such as by Larson & Whittet (2005) and 
Nishiyama et al. (2006), both of which have found a range in values for the IR exponent. 

Values for R(V) are often estimated from the formula R(V) = 1.1E(V — K)/E{B — V), which is 
based on van de Hulst's theoretical extinction curve No. 15 (e.g., Johnson 1968). In Figure [TBI we plot 
E(V — K)/E{B — V) vs. R(V) for our survey sightlines. The dotted line shows the van de Hulst relation, 
which agrees well with the data for values of R(V) near 3 — not surprising since it was derived from a 
theoretical curve with R = 3.05 — but systematically deviates at higher and lower values of R(V). The solid 
line in the figure shows the best-fit linear relation which minimizes the residuals perpendicular to the fit. It 
is given by 



EjX - V) 
E(B - V) 
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This relationship was derived only from the sightlines with 2MASS i^T-band measurements (solid circles in 
Figure I15|) but also agrees with those measurements based on Johnson Jf-band photometry (open circles) . 
This exact form of this relation depends slightly on our choice of an IR power law exponent of —1.84, and the 
small scatter is another indicator that our data are consistent with a single functional form for IR extinction. 
The relationship in Equation |8]) can be reproduced by Equation |(7j) with a wavelength of A ~ 2.1 /xm. 

6.3.3. x vs. 7 

Among the many non-correlations between extinction quantities, one of the most significant is that 
between the position of the peak of the 2175 A bump (parametrized here by xq) and its FWHM (parametrized 
by 7). The lack of a relationship between these two quantities, as first reported in Paper I, places strong 
constraints on the nature of the dust grains which produce the 2175 A feature (see, e.g., Draine 2003). 

Figure [TBI shows a plot of 7 vs. Xq for our survey sightlines, along with their l-a error bars. As in all 
previous studies, the lack of a correlation is clear. 

Figure [17] shows the distribution of bump peak positions (left panel) and widths (right panel) plotted in 
histogram form (shaded regions). With the exception of a few outliers, the distribution of bump peaks can 
be fitted well with a Gaussian function, as indicated in Figure [17] by the smooth solid curve. The centroid 
of the Gaussian is at xq = 4.5903 /im _1 and its width is given by a — 0.0191 /im -1 . These correspond to a 
mean bump position of 2178.5 A with a 1-er range of ±9.1 A. The RMS value of the xq measurement errors 
for the full sample is ±0.0058 /jm" 1 (corresponding to ±2.8 A). 

While these results suggest small but significant variations in bump positions — as reported in Paper I 
— the Gaussian-like distribution of xq values in Figure [TBJ led us to consider that perhaps our error analysis 
might underestimate the uncertainties in xq and that the width of the Gaussian itself might represent the 
true observational error. We examined this issue by considering the results for sightlines towards stars in 
open clusters and associations. We have previously used such sightlines to help estimate extinction curve 
measurement uncertainties (Massa & Fitzpatrick 1986). If it is assumed that the true extinction curve is 
identical for all cluster sightlines (based on the small spatial separation between the sightlines), then each 
cluster curve represents an independent measurement of the same curve and the variations from sightline-to- 
sightline give the net measurement errors. Since it is unlikely that there is no cosmic variability among the 
sightlines, this procedure provides upper limits on observational errors. We utilized data for the 13 clusters 
with five or more stars in our survey (NGC 457, Cr 463, NGC 869, NGC 884, NGC 1977, NGC 2244, NGC 
3293, Tr 16, NGC 4755, NGC 6231, NGC 6530, Tr 37, and Cep OB3), yielding a total of 154 sightlines. For 
each cluster we computed the mean value of xq and subtracted it from the individual cluster values. We 
then examined the ensemble of residuals for the full cluster sample. The shape of the residuals distribution is 
Gaussian-like, although very slightly skewed towards positive values. A Gaussian fit yields the result shown 
by the dotted curve in the xo panel of Figure [TBI (scaled to match the height of the main distribution), which 
has a width given by a — 0.011 /zm . This is larger than the RMS measurement error for the cluster sample, 
i.e., ±0.0058 /im _1 , but significantly smaller than the observed width of the full sample. For a number of the 
clusters, there are obvious curve-to-curve variations present, and the width of the residuals distribution must 
certainly overestimate the measurement errors. From this more detailed analysis, our conclusion remains 
that the small variations of ±9.1 A seen in the full survey sample, are significantly larger than the expected 
observational errors and indicate true variations in the position of the bump peak. 

The distribution of bump widths (right panel of Figure 1 17p is decidedly non-Gaussian, with a strong tail 
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in the direction of large values. The main peak of the distribution, however, is Gaussian in appearance and a 
fit to this region (smooth curve in the figure) yields a centroid of 0.890 /xm _1 and a width of a = 0.050 /urn - 1 . 
The RMS value of the observational errors is ±0.031 /im _1 , and so the width of this Gaussian, and the width 
of the whole distribution is clearly larger than can be accounted for by observational errors. Again consistent 
with earlier results, we find that the bump widths vary significantly from sightline to sightline, but with no 
correlation with the centroid position. The shape of the 7 distribution might suggestion two populations of 
bumps — one characterized by the Gaussian fit and the other characterized by a larger mean centroid (~1.1 
/1m -1 ) and a wider range of values (a ~ 0.1 /1m -1 ), but this is not the only possible interpretation of the 
results. We examined the bump widths for the cluster sample as above. However, the distribution of cluster 
residuals is complex, showing the asymmetry of the full sample, and indicating significant variations within 
the clusters, and we were thus not able to confirm the accuracy of the measurement errors. The range in 
observed 7 values is so large, however, that the evidence for cosmic scatter is unambiguous. 



6.4. i?(y)-Dependence 

Cardelli, Clayton, & Mathis (1988, 1989) were the first to demonstrate a link between UV and optical/IR 
extinction by showing that R(V) is related to the level of UV extinction. Essentially, sightlines with large 
R(V) values tend to have low UV extinction, and vice versa. Cardelli et al. quantified this relationship in 
the following way: 

^=a(\)+b(X)R(V)- 1 ; (9) 

i.e., the total extinction at wavelength A normalized by the total extinction at V is a linear function of 
In the time since the original work, the perception of this relationship has evolved to the point 
where it is often referred to as a "law" and Galactic extinction curves are often stated or assumed to be a 
1-parameter family (with i?(U) _1 as the parameter). Recently, for example, Valencic et al. (2004) found 
that 93% of a large sample of Galactic extinction curves obey a modified form of this relation. In this section, 
we will show that the relationship in Equation (|9|) is partially illusory and that Galactic extinction curves 
are decidedly not a 1-parameter family in i?(U) _1 . 

The original basis for Equation ^ is data such as shown in Figure [TBI where we plot A(X)/A(V) vs. 
i?(U) _1 at four different UV wavelengths for our survey sightlines (see Figure 1 in Cardelli et al. 1989). It 
is clear why a linear function would be chosen to quantify the obvious relationships seen in the figure, and 
the data give the impression of being reasonably well-correlated. The solid line in the 2695 A panel shows 
an example of such a linear relationship. It is a weighted fit which minimizes the residuals in the direction 
perpendicular to the fit, and is given by A(2695 k)/A{V) = 0.58 + 4.73 R(V)-\ The appearance of Figure 
1181 is. however, deceiving. The normalization used in the y-axis is constructed from the measured values of 
E(X - V)/E{B - V) by the transformation 

A(X) E(X-V) 

W) = E(B-V) R{V) +1 ' (10) 

Thus, the four panels in Figure [TBI essentially amount to plots of xy vs. x and, even if x and y were completely 
unrelated, some degree of apparent correlation would inevitably appear. In addition, if there actually were 
an intrinsic relationship between x and y, its significance could be greatly overinterpreted. 

The true significance of the relationships plotted in Figure [18] could be fairly assessed if the measurement 
errors were well-determined, including the effects introduced by the transformation in Equation ()10|) . Our 
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Monte Carlo error analysis allows us to quantify the uncertainties in any combination of fitted or measured 
quantities and fully take into account the artificial correlation induced by Equation (|10[) . The principal axes 
of the l-er errors for A(X) /A(V) and arc indicated in Figure IT51 by the dotted error bars and show the 

strong correlation produced by the chosen normalization. Ironically, because of the normalization and the 
resultant error correlations, the appearance of correlations in Figure[l8jis actually enhanced by uncertainties 
in R(V). When the correlated errors are taken into account, the linear relation shown in the 2695 A panel 
of Figure [TBI is found to have a reduced \ 2 value of ~4.2, indicating that the scatter about the relation is 
more than four times greater than accounted for by the known uncertainties in the data. 

Although our error analysis offers us a way to overcome the complications arising from the normalization 
chosen by Cardelli et al., a more direct approach to evaluating the relationship between UV and IR extinction 
is simply to return to the actual quantities determined by the extinction analysis and look at plots of 
E(X — V)/E(B — V) vs. R(V). A small amount of algebra (i.e., combining Equations [9] and [10] ) will show 
that, if a linear relationship as in Equation §§§ exists, then the following relation must also hold: 

lf B ~_ V v) =b(\) + [a(X)-l] R(V) . (11) 

I.e., any true linear correlation between A(X)/A(V) andi?(V) _1 will also be present between E(X—V)/E(B — 
V) and R(V), with a simple transformation relating the linear coefficients. The great advantage to viewing 
the data in the E(X — V)/E(B — V) vs. R(V) reference frame is that no artificial correlations (in either 
the parameters or their uncertainties) are introduced and we only have to deal with the natural correlations 
which arise from the dependence of all the extinction parameters on the properties of the best-fit stellar 
SEDs. Figure H9l illustrates this point. It shows plots of E(X - V)/E(B - V) vs. R(V) (along with their 
associated 1-a errors) for the same four wavelengths as in Figure [IB] The 2695 A panel shows a weighted 
linear fit, which is given by £'(2695 k-V) / E(B - V) = 4.80 - 0.44 R(V). This is nearly exactly what would 
be expected from the fit to A(2695 A)/A(V) vs. i?(V) _1 in Figure [TBI and the transformation in Equation 
(llip . The exact transformation is shown by the nearly coincident dashed line in Figure [T9l verifying the 
argument leading to Equation (jTTJl . (A careful comparison between the lines and the data points will show 
that the two panels do indeed show the same relationship.) 

From a statistical point of view, the presentations in Figures [TB] and [TH] are identical, as are the linear 
fits to the 2965 A data. In fact, both fits have nearly identical, underwhelming, reduced \ 2 values of ~4.2. 
Fits performed at wavelengths shortward of 2675 A, show increasingly large values of \ 2 (e-g-, X 2 = 6.9 for 
a linear fit at 1665 A). The lesson from Figures [TBI and [T9l is that the choice of normalization affects the 
perception of how well-related R(V) is to the level of UV extinction. The data in Figure [15] look better- 
correlated than do those in Figure [TH — but they are not. The "cosmic scatter" is appreciable and there 
is no functional relationship between R{V) and the UV extinction (whether linear or more complex) for 
which the scatter approaches the current level of measurement errors. Although there is a relationship in the 
sense that large- R(V) curves [i.e., R(V) > 4.0] differ systematically from low-i?(V) curves, UV extinction 
properties cannot be expressed as a 1-paramctcr family in R(V) at anywhere near the level of observational 
accuracy. 

Note that the discussion above is not affected by the likelihood that some of the sightlines are composites, 
possibly spanning distinct regions of very different R(V). If a universal relationship of the form in Equation 
© were to hold, then composite regions would still lie along the line a(A) -I- 6(A) i?(V) _1 and would not 
result in "cosmic scatter." 



The recent study by Valencic et al. (2004) examined the i?(V)-dependence of UV extinction by looking 
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at correlations between i?(V) _1 and UV fit parameters based on the Paper III formulation, as had been done 
in F99 and F04. As noted earlier, they concluded that most Galactic curves (i.e., 93%) are consistent with a 
i?(V)-dependent "law." While the approach was somewhat different, the results of the Valcncic ct al. study 
suffer the same problem as shown in Figures IT51 and [TO1 because they explicitly multiply the UV fit coefficients 
by i?(V A ) _1 , thus forcing a correlation in the errors and enhancing the perception of a correlation between 
the quantities. The large percentage of curves believed to be consistent with a single i?(U)^ 1 relation results 
from the complications introduced by the choice of curve normalization. Figure [20] shows the relationship 
between R(V) and the UV fitting parameters in their original form, in which the correlations among the 
parameters and their errors are not artificially enhanced. The figure contains several extra quantities derived 
from the fit parameters which help describe features of the UV curves and which are defined in Sj3] (and also 
in the figure legend). This figure again shows that, while systematic differences exist between high-i?(V) 
and low-i?(V) curves, there is no simple (or complex) relation between R(V) and the UV fitting properties 
which is consistent with measurement errors. 



7. SUMMARY 

Several of our findings are worth emphasizing. 

Variability of xq and 7: We found, in accordance with our previous analysis, that while the central 
position of the 2175 A bump does not vary much, it is indeed variable. On the other hand, the bump width 
varies considerably. Further, there still appears to be no relationship between the central position of the 
bump and and its width - verifying the results described in Paper I to a higher degree of accuracy than 
previously possible. 

Correlations among the fit parameters: Generally, the various fit parameters which describe the shape 
of the UV-through-IR extinction curves are not related to one another. The only exceptions are c\ — c-i and 
kiR — R(y)- The first relates the slope and intercept of the UV portion of the curve and is effectively 
a functional relationship for most of the data. The latter relates the scale factor of the power law used 
to describe the IR portion of the curve and the ratio of total-to-selective extinction. At the level of our 
measurements errors, these two parameters are functionally related and consistent with a universal form to 
IR extinction. However, due to observational uncertainties, the dominance of diffuse ISM sightlines in our 
study, and the limitations of the IR data we employ, our results are not ideal for addressing the detailed 
shape of IR extinction and significant sightline-to-sightline variations could still exist. 

Correlations with dust density: As determined in the past, various properties of the extinction curves 
are weakly correlated with the mean line of sight dust density, as measured by E(B — V)/d. Presumably, 
these correlations reflect the operation of a physical process, such as grain growth or coagulation. 

Correlations of R(V) with UV curve properties relationship: Correlations between R(V) and UV 
curve properties have received considerable attention over the years since Cardelli et al. (1988) first pointed 
out that a large-scale trend could be found. By expressing the curves in their native form, we verify that 
a there is a weak relationship between R(V) and the UV in the sense that sightlines with extremely large 
R(V) values tend to have low normalized UV extinction curves. However, this relationship is only evident 
for the largest R(V) sightlines. For the majority of sightlines, there is no evidence for an i?(V)-dependence 
of the extinction curve shapes. Specifically, inspection of Figure Q1J] shows that for diffuse ISM sightlines 
(2.4 < R(V) < 3.6), which comprise the bulk (82%) of our sample, no relationship exists, even though a 
large range in the extinction is present. Moreover, Figure l20l shows that the 2175 A bump properties display 
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a similar behavior, i.e., a trend is evident only for the largest R(V) values, with the strength of the bump 
(as measured by Ab ump or Eb ump ) tending to be slightly weaker than the average for lower R(V) sightlines. 
Thus, we conclude that there is no global 1-parameter family of extinction curves, although extremely large 
RiV) curves tend to have distinctive properties. 

Extinction curve variability the meaning and utility of an average curve: The previous result 
begs the meaning of an average extinction curve. As noted in ij6.U simple mean curves always reflect the 
biases of their parent samples. Our mean curve for the diffuse ISM (Figure [HI an average of all curves with 
2.4 < R(V) < 3.6, and that from Valencic et al. (2004, for R=3.1) are derived from the largest samples and 
probably provide the best estimate of mean Galactic diffuse ISM extinction properties at short wavelengths. 
However, one must always be mindful of the dark shaded region in Figure [5] which illustrates the RMS 
variance that can be expected for an extinction curve along an arbitrary diffuse sightline. Typical RMS 
dispersions in fc(A - V) are 0.31, 0.68, 0.62, 1.44 at A =2695, 2175, 1665, 1245 A, respectively. This means 
that if the mean curve is used to deredden an object with E(B — V) = 0.50 mag, the uncertainty in the 
dcrcddened continuum at these wavelengths would be 0.15, 0.30, 0.30 and 0.72 mag, respectively, due to 
uncertainties in the extinction alone! It is, however, possible to take advantage of localized uniformity in the 
extinction to reduce this error (e.g., Massa & Savage 1981). 

Physical implications: Perhaps the best way to summarize the physical origin of the i?(V)-dependence 
for Galactic extinction is to examine the i?(V^)~ 1 — C2 plot shown in Figure [2T1 which includes a best-fit 
linear relation (such a relation was the basis for the i?(F)-dependent curves produced by F99 and F04 - see 
Figure 10 of F04). Figure [21] can be summarized thusly: when extinction curves are steep in the optical 
(large i?(V^) _1 ) they tend to stay steep in the UV (large cq) and when extinction curves are flat in the 
optical (small i?(U) _1 ) they tend to stay flat in the UV (small o%). However, as the scatter in the Figure 
illustrates, this trend is only apparent for extreme values of R{V}~ X or C2- In general, there is no unique 
relation between these parameters over the range spanned by most of the sample. It is likely that the general 
connection between UV and optical extinction slopes simply reflects the fact that the overall grain size 
distribution affects all wavelengths. But the presence of such a large scatter demonstrates that several other 
factors (e.g., chemical composition, grain history, coagulation, coating, radiation environment, etc.) must 
also be involved. In other words, the large variance in the relation between UV and optical slopes indicates 
that dust grain size distributions do not behave as a 1-parameter family. 

Final remarks: Having painted a negative picture of the relationship between R(V) and extinction prop- 
erties at other wavelengths, it would be disingenuous of us to present yet another set of i?(F)-dependent 
curves. The results of Cardelli et al. (1989), F04 (which supersedes those of F99), and Valencic et al. 
(2004) are all reasonable, and the differences among them are instructive of the biases introduced by sam- 
ple selection and methodology. It should always be remembered that these curves represent very general 
trends in Galactic extinction and do not constitute a standard of normalcy. Finally, although the prospect 
of accurately dereddening an unknown SED using a mean extinction curve is poor, this same variability 
demonstrates that extinction curves are responsive to local conditions, so that each one contains potentially 
unique information about the grains along the sightline. 
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A. Obtaining the Data from this Study 

The extinction curves and associated information from this study are available via anonymous ftp at 
"ftp.astronomy.villanova.edu". After logging in as "anonymous", change to the appropriate directory by 
typing "cd fitz/FMV.EXTINCTION". A "README" file, two IDL procedure files, and two tar files are 
present. The tar files contain the extinction curves from this program and the R-dependent curves from 
Fitzpatrick 2004 ("FMV.EXTCURVES.tar" and "F04_RCURVES.tar" , respectively). The tar files unfold 
into individual ascii data files for each extinction curve, containing the 12 extinction parameters which 
describe the curve and the curve itself (along with 1-a uncertainties for the curves from this paper). The 
two IDL procedures provide tools for reading data from the ascii files and for constructing extinction curves 
from the 12 fit parameters ("READ_FMV_FILES.pro" and "MAKE_FMV_CURVE.pro" , respectively). 
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Table 1 . Basic Data for Survey Stars [The complete version of this table is in the electronic edition of the 

Journal. The printed edition contains only a sample.] 



Star a 


Spectral 


V 


Distance 


I 


6 


Reference 




Type" 


(mag) 


(PC) 


(°) 


(°) 




HD698 


B5 II: SB 


7.10 


1125 


117.689 


-4.25 


1 


HD3191 


Bl IV:nn 


8.58 


1203 


121.068 


-1.36 


2 


BD+57 245 (NGC 457) 




9.85 


2429 


126.583 


-4.58 




BD+57 252 (NGC 457) 


Bl IV 


9.51 


2429 


126.644 


-4.42 


2 


NGC 457 Pcsch 34 




10.61 


2429 


126.646 


-4.38 




NGC 457 Pcsch 13 




10.78 


2429 


126.646 


-4.38 




NGC 457 Pcsch 9 


Bl V 


9.83 


2429 


126.646 


-4.38 


3 


Cr 463 #18 




10.35 


702 


127.091 


9.20 




BD+70 131 (Cr 463) 




10.06 


702 


127.280 


9.17 




Cr 463 #5 




10.37 


702 


127.264 


9.37 





a The stars arc listed in order of increasing Right Ascension using the most commonly adopted 
forms of their names. The first preference was "HDnnn" , followed by "BDnnn" , etc. There are 
185 survey stars which are members of open clusters or associations, or 56% of the sample. The 
identity of the cluster or association is either contained in the star name itself (e.g., NGC 457 
Pesch 34) or is given in parentheses after the star's name. 

b Spectral types were selected from those given in the SIMBAD database, and the source of 
the adopted types is shown in the "Reference" column. When multiple types were available for a 
particular star, we selected one based on our own preferred ranking of the sources. 

c The NGC 2244 distance is from Perez et al. 1987; the NGC 3293 distance is from Bolona & 
Crampton 1974; the Trumplcr 14 and 16 distances are from Masscy & Johnson 1993; the Ccp OB3 
distance is from Crawford & Barnes 1970. The distances to all other clusters or associations are 
from the Open Clusters and Galactic Structure database maintained by Wilton S. Dias, Jacques 
Lcpine, Bruno S. Alcssi, and Andre Moitinho at http://www. astro. iag.usp.br~wilton/ For the 
non-cluster stars, distances were calculated using the E(B — V) values from this study and the 
absolute magnitudes from Turner 1980 (for mid-B and earlier types) and Blaauw 1963 (for mid-B 
and later types). 

References. — (1) Hiltner 1956; (2) Morgan, Code, & Whitford 1955; (3) Hoag, & Applequist 
1965; (4) Mai'z-Apellaniz, Walborn, Galue, & Wei 2005; (5) Johnson, & Morgan 1955; (6) Slettebak 
1968; (7) Schild 1965; (8) Racine 1968; (9) Morgan, Keenan, & Kellman 1943; (10) Mendoza 1956; 
(11) Morgan, Whitford, & Code 1953; (12) Osawa 1959; (13) Cowley, Cowley, Jaschek, & Jaschek 
1969; (14) Metreveli 1968; (15) Guetter 1968; (16) Boulon, & Fehrenbach 1959; (17) Buscombe 
1962; (18) Roman 1955; (19) Georgelin, Georgelin, & Roux 1973; (20) Bouiguc 1959; (21) Wenzel 
1951; (22) Smith 1972; (23) Penston, Hunter, & O'Neill 1975; (24) McNamara 1976; (25) Levato 
H., & Abt H.A.; (26) Schild, & Chaffee 1971; (27) Borgman 1960; (28) Warren, & Hesser 1977; (29) 
Sharplcss 1952; (30) Racine 1968; (31) Crawford, Limber, Mendoza, Schultc, Stcinman, & Swihart 
1955; (32) Johnson, & Morgan 1953; (33) Meadows 1961; (34) Johnson 1962; (35) Hoag, & Smith 
1959; (36) Barbier, & Boulon 1960; (37) Young 1978; (38) Claria 1974; (39) Feast, Thackeray, & 
Wesselink 1955; (40) Moffat, & Fitzgerald 1974; (41) Houk 1978; (42) Feast, Stoy, Thackeray.fe 
Wesselink 1961; (43) Denoyellc 1977; (44) Hoffleit 1956; (45) Turner, Grieve, Herbst, & Harris 
1980; (46) Morrell, Garcia, & Levato 1988; (47) Masscy, & Johnson 1993; (48) Levato, Malaroda, 
Morrell, Garcia, & Hernandez 1991; (49) Morris 1961; (50) Hiltner, Garrison, & Schild 1969; (51) 
Houk, & Cowley 1975; (52) Seidensticker 1989; (53) Feast, Thackeray, & Wesselink 1957; (54) 
Perry, Franklin, Landolt, & Crawford 1976; (55) Schild 1970; (56) Do Vaucouleurs 1957; (57) 
Walraven, & Walraven 1960; (58) Garrison 1967; (59) Buscombe 1969; (60) Hardie, & Crawford 
1961; (61) Perry, Hill, & Christodoulou 1991; (62) Schild, Ncugebaucr, & Westphal 1971; (63) 
Schild, Hiltner, & Sanduleak 1969; (64) Levato, & Malaroda 1980; (65) Garrison, & Schild 1979; 
(66) Roman 1956; (67) Houk, & Smith-Moorc 1988; (68) Hiltner, Morgan, & Ncff 1965; (69) 
Walker 1957; (70) Walker 1961; (71) Stcbbins, & Kron 1956; (72) Guetter 1964; (73) Hill 1970; 
(74) Hcrbig, & Spalding 1955; (75) Vrba, & Rydgrcn 1984; (76) Hack 1953; (77) Roman 1951; 
(78) Divan 1954; (79) Garrison, & Kormendy 1976; (80) Simonson 1968; (81) Blaauw, Hiltner, & 
Johnson 1959; (82) Pecker 1953; (83) Boulon, Duflot, & Fehrenbach 1958 
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Table 2. Adopted Temperature Scale for Main Sequence O Stars 



Spectral 


T eff 


Type 


(K) 


06 


40000 


06.5 


38500 


07 


37000 


07.5 


36500 


08 


36000 


08.5 


34750 


09 


33500 


09.5 


32750 


BO 


32000 
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Table 3. Best-Fit Parameters for Survey Stars [The complete version of this table is in the electronic 
edition of the Journal. The printed edition contains only a sample.] 



Star 




logg b 


[m/H] c 


^turb 


Br 


E(B-V) 




(K) 






(km/s) 


(mas) 


(mag) 


HD698 


18434 ± 499 


3.72 ±0.31 


-0.18 ±0.10 


3.2 ± 0.5 


0.0670 ± 0.0014 


0.37 ±0.01 


HD3191 


24001 ± 890 


3.41 ±0.28 


-0.43 ±0.06 


10 


0.0328 ± 0.0009 


0.68 ±0.01 


BD+57 245 (NGC 457) 


22885 ± 697 


3.64 ±0.20 


-0.39 ±0.08 


5.5 ± 0.6 


0.0160 ± 0.0004 


0.50 ±0.01 


BD+57 252 (NGC 457) 


24924 ± 616 


3.64 ±0.16 


-0.47 ±0.06 


10 


0.0173 ± 0.0004 


0.51 ±0.01 


NGC 457 Pesch 34 


23594 ±661 


3.91 ±0.17 


-0.07 ±0.06 


1.2 ± 0.8 


0.0110 ±0.0003 


0.51 ±0.01 


NGC 457 Pesch 13 


22023 ± 809 


3.84 ±0.17 


-0.47 ±0.11 


2.1 ± 0.8 


0.0114 ±0.0003 


0.51 ±0.01 


NGC 457 Pesch 9 


25738 ± 577 


3.77 ±0.16 


-0.75 ±0.07 


9.1 ± 0.7 


0.0147 ± 0.0003 


0.54 ±0.01 


Cr 463 #18 


11891 ± 286 


3.99 ±0.16 


-0.60 ±0.10 





0.0198 ± 0.0003 


0.35 ±0.01 


BD+70 131 (Cr 463) 


11351 ± 204 


3.93 ±0.17 


-0.73 ±0.07 


5.1 ± 0.4 


0.0210 ± 0.0003 


0.28 ±0.01 


Cr 463 #5 


11859 ± 230 


4.10 ±0.19 


-0.59 ±0.09 


2.0 ±0.4 


0.0170 ± 0.0002 


0.30 ±0.01 



a For the O stars analyzed using the TLUSTY atmosphere models, the values of T e f f were adopted from the spectral 
type vs. T e ff relation given in Table 2. These stars can be identified by their l-cr uncertainties, which are ±1000 K. 

b For stars in clusters, the surface gravities are determined as discussed in §3.1 and rely on stellar evolution models 
and cluster distance determinations. Surface gravities for non-cluster stars are not always well-determined, because 
of a lack of specific spectroscopic indicators. In some cases, the best-fit solutions for these stars indicated physically 
unlikely results (i.e., logg > ~4.3 or logg < ~3.0). For these stars, a value of logg = 3.9 was assumed (which is the 
mean logg of the rest of the sample) and a 1-<t uncertainty of ±0.2 was incorporated in the error analysis. These 
cases can be identified by logg entries of "3.9 ± 0.2." 

c For the O stars in the sample, our fitting procedure utilized solar abundance TLUSTY models. For these stars the 
values of [m/H] are indicated by entries of "0" without uncertainties. 

d For the O stars, the adopted TLUSTY models incorporate v tur i, = 10 kms -1 . For these stars the values of Vturb 
are indicated by entries of "10" without uncertainties. For the B stars, which were modeled using ATLAS9 models, 
the values of Vf ur i were determined by the fitting procedure, but were constrained to lie between and 10 kms -1 . 
Stars whose best-fit SED models required these limiting values are indicated by u tIlr f, entries of "0" or "10" , without 
error bars. The uncertainties for stars with best-fit Vt ur b values close to these limits may be underestimated due to 
this truncation. 



Table 4. Best-Fit Extinction Curve Parameters for Survey Stars [The complete version of this table is in the electronic edition of the 

Journal. The printed edition contains only a sample.] 











UV Coefficients 








Optical Spline Points b 


IR Coefficients 


Star 


x 


7 


Cl 


C2 


C3 


C 4 


C5 


Oi o 2 o 3 


R(V) klR 



HD698 


4, 


551 


± 





.006 





.96 


± 


0. 


.03 





.07 


± 





19 





.99 


± 





05 


2 


95 


± 





.19 


0. 


15 ± 





.05 





.51 ± 





.29 


2.38 


± 


0.09 


1 


.33 


0. 


.00 


3 


.94 


± 





.16 


1. 


70 ± 0. 


18 


HD3191 


4. 


636 


± 





.003 





.94 


± 





.02 


-0 


.79 


± 





20 


1 


.00 


± 





04 


2. 


99 


± 





16 





31 ± 





.02 


5 


.76 ± 





.09 


2.04 


± 


0.06 


1 


.33 


0. 


.01 


2 


.81 


± 





.09 


0. 


93 ± 0. 


14 


BD+57 245 (NGC 457) 


4. 


561 


± 





.005 





.89 


± 


0. 


.03 


-0 


.47 


± 


0. 


.31 





.88 


± 


0. 


06 


3. 


.13 


± 





.19 


0. 


.24 ± 





.03 


5 


,71 ± 





.16 


2.05 


± 


0.06 


1 


.31 


0. 


00 


2 


.97 


± 


0. 


13 


0. 


.96 ± 0. 


17 


BD+57 252 (NGC 457) 


4. 


577 


± 





.005 





.92 


± 


0. 


.02 


-0 


.66 


± 





.25 





.95 


± 





.05 


3 


39 


± 





.17 


0. 


17 ± 





.02 


5 


.48 ± 





.23 


2.15 


± 


0.06 


1 


.32 


0. 


.00 


2 


.97 


± 


0. 


.10 


0. 


98 ± 0. 


15 


NGC 457 Pesch 34 


4. 


579 


± 





.003 





.91 


± 


0. 


.01 


-0 


.55 


± 


0. 


.20 





.86 


± 





04 


3 


32 


± 





11 


0. 


31 ± 





.02 


6 


.14 ± 





.09 


1.99 


± 


0.07 


1 


.30 


0. 


01 


2 


.94 


± 


0. 


.11 


0. 


.84 ± 0. 


16 


NGC 457 Pesch 13 


4. 


587 


± 





.004 





.87 


± 


0. 


.02 


-0 


.09 


± 


0. 


33 





.84 


± 


0. 


.07 


2. 


91 


± 





.14 


0. 


26 ± 





.02 


5 


.59 ± 





.17 


2.03 


± 


0.07 


1 


.31 





01 


3 


.11 


± 


0. 


.11 


0. 


90 ± 0. 


16 


NGC 457 Pesch 9 


4. 


578 


± 





.006 


1 


.07 


± 


0. 


03 


-1 


.13 


± 


0. 


22 


1 


.08 


± 


0. 


.04 


4. 


.31 


± 





.22 





18 ± 





.02 


5 


.00 ± 





.18 


2.37 


± 


0.05 


1 


.35 





.01 


2 


.76 


± 


0. 


.07 





.79 ± 0. 


12 


Cr 463 #18 


4. 


606 


± 





.009 


1 


.18 


± 





.05 


-0 


.47 


± 





26 


1 


.03 


± 


0. 


.07 


5. 


23 


± 





.45 





30 ± 





.03 


5 


,09 ± 





.20 








1 


.37 


0. 


01 


3 


.38 


± 





.17 




1.29 




BD+70 131 (Cr 463) 


4. 


577 


± 





.011 


1 


11 


± 





.05 





.17 


± 


0. 


.32 





.86 


± 





08 


5. 


11 


± 





.49 





35 ± 





.04 


5 


.00 ± 





19 








1 


.37 


0. 


00 


3 


.40 


± 





20 




1.30 




Cr 463 #5 


4, 


610 


± 





.008 


1 


.08 


± 


0. 


.04 





.45 


± 


0. 


33 





.74 


± 


0. 


08 


4. 


.77 


± 





.34 


0. 


.32 ± 





.03 


5 


.06 ± 





.18 








1 


.37 


0. 


00 


2 


.91 


± 


0. 


.18 




0.99 





a For the stars HD237019, HD18352, and HD25443 the long wavelength IUE spectra are incomplete. For these cases we constrained the UV linear extinction component to follow the relation 
ci — 2.18 — 2.91C2 from Fitzpatrick 2004. For these stars we list uncertainties for the C2 values but not for the c± values. 

b The uncertainties in the O2 and O3 optical spline points (at wavelengths of 4000 A and 5530 A, respectively) are typically 0.01 or less and arc not listed. For several stars — those without U band 
photometry — we did not solve for the 0± point at 3300 A. , 

c For field stars without IR photometry, we assumed ^CV") — 3.1 and km — 1.11, with the latter based on the relation km. — 0.63i?(V) — 0.84 from Fitzpatrick 2004. For such st^s in clusters, we 
adopted the mean R(V) of the other cluster members and a value of km based on the aforementioned relation. These assumed values are listed in the Table without uncertainties. Several survey 
stars have apparently noisy JHK data and yielded very uncertain values of km. For these, we ultimately derived the extinction curve by solving for the best-fit value of R(V) with kjR constrained to 
follow the Fitzpatrick 2004 relation. The resultant R(V) values are listed with their uncertainties while the k IR values are listed without uncertainties. 



-33 - 



Table 5. Extinction Curve Parameters for Mean Curve in Figure 9 



Parameter 


Value 


x 


4.592 /im" 1 


7 


0.922 ^nT 1 


Cl 


-0.175 


C2 


0.807 


C3 


2.991 


C 4 


0.319 


c 5 


6.097 


Oi 


2.055 


o 2 


1.322 


o 3 


0.000 




3.001 


km 


1.057 



+ 90 



90' 



Cr463 / Tr37 
C'epOB3 



o«-#- 



t h & x NGC457 



NGC4755 

Carina 



Fig. 1. — Sky distribution of survey stars, shown in a sinusoidal projection of Galactic coordinates. The 
Galactic plane is the solid horizontal line and the Galactic center is in the middle of the figure. The locations 
of open clusters or associations containing five or more survey stars are indicated by the large circles, labeled 
with the cluster name. The sizes of the circles do not represent the physical extent of the clusters or 
associations. 
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Fig. 2. — Representative properties of the survey stars. The three panels show results from the analysis 
presented in §4l which help characterize the properties of our sample. The typical survey member is a 
micl-B star (top panel), with a median reddening of E(B — V) = 0.45 mag (middle panel), viewed along a 
line-of-sight passing through the diffuse ISM (bottom panel). The bottom panel also shows a Gaussian fit 
to the values of R(V) = Ay / 'E(B-V) in the neighborhood of the peak in the distribution. The peak is 
located at R{V) — 2.99, similar to the mean values usually attributed to the diffuse ISM, and the width of 
the Gaussian corresponds to a = ±0.27. The mean and median values of R{V) for the whole sample are 
3.22 and 3.05, respectively. 
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Fig. 3. — A parametrized representation of normalized UV-through-IR extinction (solid and dashed curve). 
The curve consists of two parts: 1) A < 2700 A (shaded region), where we adopt a modified version of 
the 3-component parametrization scheme of Paper III; and 2) A > 2700 A, where we adopt a cubic spline 
interpolation through sets of IR (I n ), optical (O n ), and UV (U n ) "anchor points." The values of the anchor 
points and the seven parameters describing the UV portion of the curve are determined by fitting the observed 
SED of a reddened star, as described in |J3J The particular curve shown in this figure corresponds to that 
derived for the star HD147933. 
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Fig. 4. — Normalized extinction curves for 328 Galactic stars derived using the extinction-without-standards 
approach. The symbols show the normalized ratios of the model atmosphere fluxes to 1) the IUE spec- 
trophotometry in the UV (A -1 > 3.3 /im -1 ), 2) Johnson UBV photometry in the optical, and 3) Johnson 
or 2MASS JHK photometry in the near-IR (A -1 < 1 /im -1 ). Individual l-a error bars are shown for the 
data points, but are typically only visible in the region of the 2175 A bump and in the far-UV, where the 
signal level of the IUE data is lowest. Small crosses indicate IUE data points excluded from the fit, for 
the reasons discussed in Paper IV. IUE data points in the region 1215 < A < 1235 A are excluded due to 
contamination from scattered solar Lya photons. The solid curves are the parametrized fits to the data as 
determined by the SED-fitting procedure discussed in <J31 and the shaded regions show the l-a uncertainty 
in the curves, based on Monte Carlo simulations. For comparison, the dash-dot curves show an estimate of 
the average Galactic extinction curve from F99 (corresponding to R(V) — 3.1). Only the first panel of the 
figure is shown in the print edition of the Journal. The entire figure consists of 33 panels and is available in 
the electronic edition. 
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Fig. 5. — logTeff as a function of spectral type for the survey stars of type BO or later. Luminosity class 
V or IV stars are shown by rilled circles; higher luminosity classes by open circles. For comparison, data for 
the 44 unreddened stars from the photometric calibration study of Fitzpatrick & Massa (2005b) are shown 
by the open squares. The various dotted, dashed, and solid lines show a number of published spectral type 
vs. T e ff calibrations, as indicated in the figure. Small random horizontal offsets have been added to the 
data points to increase their visibility. 
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Fig. 6. — UV spectra of two survey stars, HD228969 and BD+54 973, compared with three early-type 
spectral classification standard stars. The reddened survey stars are classified as types B2 II: and B3 V, 
respectively, although the SED-fitting procedure indicates temperatures in the neighborhood of 30000 K. 
Comparison of the UV spectral features with the classification standards indicates that these two stars are 
likely misclassified members of earlier spectral classes. Much of the scatter in the T e g vs. spectral type 
diagram in Figure [5] is probably the result of such classification uncertainties. 
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Fig. 7. — T e g as a function of the Stromgren reddening-free index [c] = c\ — 0.20(6 — y), which measures the 
strength of the Balmer jump, for 162 sample stars. Symbols are the same as for Figure [5j with the addition 
of the open circles which denote O-type stars. Also shown are the spectral type vs. T c g calibrations from 
Figure (3 transformed into the T e g vs. [c] plane using Crawford's (1978) spectral type vs. cq and cq vs. 
(b — y)o relations. Note that Stromgren photometry was not used in modeling the SEDs of our survey stars. 
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Fig. 8. — Analytical representations of 328 Galactic UV-through-IR extinction curves plotted in their native 
normalization E(X — V)/E(B — V) (top panel) and transformed into an IR-normalized form A(X)/A(l/im) 
(bottom panel). The curves, in their standard form, can be reproduced from the parameters given in Table 
|U The convergence of the curves in the optical region (top panel) and the IR region (bottom panel) results 
from their normalizations. The value of R(V) = A(V)/E{B — V) for each curve is the negative of its intercept 
at A" 1 = ^m" 1 in the E{\ - V)/E(B ~ V) panel. 
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Fig. 9. — An average UV-through-IR extinction curve for our sample compared with other mean extinction 
curves. The thick solid curve is the mean for the 243 stars in our sample with 2.4 < R(V) < 3.6, i.e., for those 
sightlines with R(V) values indicative of the diffuse ISM. The dark gray shaded region shows the sample 
variance about the mean curve at all wavelengths. The light grey shaded area shows the larger variance that 
results when we include all 298 sightlines with measured R(V) values. The mean curve can be reconstructed 
from the extinction parameters listed in Table [5] The dashed and dotted curves show estimates of mean 
UV Galactic curves from the sources indicated in the figure. The large filled circles are means from the 
ANS satellite extinction catalog of Savage et al. (1985) for 800 stars with E(B — V) > 0.20 mag. These 
measurements result from filter photometry centered at wavelengths of 3000, 2500, 2200, 1800, and 1550 A. 
Error bars on the ANS data show the sample variances and include the effects of spectral mismatch. 
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Fig. 10. — Unweighted mean UV extinction curve for 318 survey stars (open circles). Ten stars (HD14092, 
BD+56 501, HD14321, HD18352, HD25443, CD-42 4819, HD99872, HD326328, HD197702, and HD239710) 
were excluded from the mean because they have incomplete IUE spectra. The individual data points (shown 
in Figured]) are typically 10 A apart. The solid curve is a parametrized fit to the mean curve, using the 
formulation discussed in Sj3j The O — C residuals are shown as filled points; they are offset from their zero 
mean for display. A number of distinct features are seen in the residuals and labeled in the figure. These 
arise from ISM gas-phase absorption lines, mismatch between the C iv A1550 stellar wind lines in the O 
stars and the static model SEDs, and known inadequacies in the ATLAS9 opacity distribution functions 
(labeled "b" in the figure; see Paper IV). No other credible features are seen in the residuals. The standard 
deviation of the residuals about their mean value of zero (excluding the labeled points) is 0.06E(B — V) mag, 
corresponding to ~ 0.02^4(V) mag. This figure demonstrates both the intrinsic smoothness of UV extinction 
and the ability of our parametrization scheme to reproduce the shape of UV extinction curve to extremely 
high accuracy. 
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Fig. 11. — Spatial trends in extinction properties. The values of four parameters which describe the shapes 
of interstellar extinction curves are plotted against Galactic longitude for each of our survey sightlines. Small 
filled circles show field sightlines. Other symbols denote sightlines to clusters or associations for which five 
or more members are included in our survey. From left to right, the clusters are: h & \ Per (x's); NGC 457 
(open circles); Cr 463 (filled diamonds); Cep OB3 (filled squares); Trumpler 37 (filled triangles); NGC 6530 
(open squares); NGC 6231 (*'s); NGC 4755 (open diamonds); Carina clusters (large filled circles, includes 
Trumpler 14 and 16, Cr 228, and NGC 3293); NGC 1977 (open triangles); and NGC2244 (+'s). The dashed 
lines show the parameter values corresponding to the diffuse mean ISM curve from Figure [Sj 
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Fig. 12. — Density trends in extinction properties. The values of four parameters which describe the shapes of 
interstellar extinction curves are plotted against E(B — V) /d for our survey sightlines. Small filled circles are 
field sightlines. Other symbols (same as in Figure fTTj) denote sightlines to clusters or associations containing 
five or more survey stars. The dashed lines show the parameter values corresponding to the diffuse mean 
ISM curve in Figure M M 



-46 - 




0.0 0.5 1.0 1.5 2.0 

c2 (UV slope) 



Fig. 13. — Plot of the slope c 2 vs. the intercept c\ for the linear component of UV extinction. The 
uncertainties in c\ and C2 are strongly correlated, as indicated by the \-a error bars (dotted lines). The 
orientation of the error bars was determined by fitting ellipses to the loci of results from the Monte Carlo 
error simulations. The individual sets of error bars are not orthogonal to each other because the scales of 
the x- and y-axes are different. The solid line is a weighted linear fit to the data, designed to minimize the 
residuals in the direction perpendicular to the fit. It is given by c\ = 2.09 — 2.84c 2 . 
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Fig. 14. — Plot of R(V) vs. the IR scale factor km (see Eq. [4]). The uncertainties in R(V) and km are 
strongly correlated, as indicated by the 1-a error bars (dotted lines). The orientation of the error bars was 
determined by fitting ellipses to the loci of results from the Monte Carlo error simulations. The solid line is 
a weighted linear fit to the data, designed to minimize the residuals in the direction perpendicular to the fit. 
It is given by k IR = -0.83 + 0.63i2(V). 
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Fig. 15.— Plot of E(V - K)/E{B - V) vs. R{V) for sightlines with 2MASS if-band photometry (filled 
circles) and Johnson if-band photometry (open circles). The dotted line shows the widely-used relation 
R{V) = 1.1E(V - K)/E{B - V). The solid line is a linear fit to the 2MASS data, minimizing the residuals 
in the direction perpendicular to the fit. It is given by R(V) = -0.26 + 1.19E(V - K)/E{B - V). 
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Fig. 16. — Plot of the FWHM of the 2175 A bump 7 vs. its central position, x , both in units of ^m _1 . 
Dotted lines show the l-a error bars, revealing the uncertainties in 7 and xq to be generally uncorrelated. 
The orientation of the error bars was determined by fitting ellipses to the loci of results from the Monte 
Carlo error simulations. The individual sets of error bars are not orthogonal to each other because the scales 
of the x- and y-axes are different. A number of sightlines which deviate most from the main distribution 
have been labeled. 
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Fig. 17. — 2175 A bump statistics. The left panel shows a histogram of the distribution of the 2175 A 
bump peak positions for our sample. The smooth curve is a Gaussian fit centered at xo — 4.5903 jUm 
(2178.5 A) with a width of a = 0.0191 /im _1 (9.1 A). The RMS value of the measurement errors in x$ is 
±0.0058 fim^ 1 (±2.8 A). The dashed curve shows a Gaussian fit to the distribution of bump positions for 
154 stars located in the 13 open clusters with more than 5 stars in the survey. For each cluster, a mean 
value of Xg was computed and the distributions constructed relative to the mean (see H6.3.3[) . The width of 
the cluster Gaussian is given by a — 0.011 fim^ 1 (5.2 A). The right panel shows the distribution of bump 
FWHM values. The smooth curve is a Gaussian fit to the main peak in the distribution. Its central position 
corresponds to 7 = 0.890 /im _1 and its width is a = 0.050 /zm . The RMS measurement error in 7 is 
±0.031 pxar 1 . 
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Fig. 18. — Plots of the normalized extinction values A{\) / A(V) vs. Rty)^ 1 at four different UV wavelengths. 
A(X) / A(V) is derived from the measured values of E(X — V)/E(B — V) by using Equation (| 10[) . 1-a error 
bars are shown, based on our Monte Carlo simulations. The error bars, which are not orthogonal to each 
other because the scales of the x- and y-axes are different, show the clear correlation in the uncertainties for 
A(X)/A(V) and R(V)-\ The straight line in the 2695 A panel is a weighted linear fit to the data, which 
minimizes the scatter perpendicular to the relation. It is given by ,4(2695 k)/A(V) = 0.58 + 4.73 R(V)-\ 
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Fig. 19. — Plots of the measured extinction values E(X — V)/E(B — V) vs. R(V) at four different UV 
wavelengths. 1-cr error bars are shown, based on our Monte Carlo simulations as in previous figures. Any 
true relationships seen in plots of A(X)/A(V) vs. are preserved in plots of E(X — V)/E(B — V) vs. 

R(V), but with the great benefit that the correlations between the errors are greatly reduced. The solid line 
in the 2695 A panel is a weighted linear fit to the data, which minimizes the scatter perpendicular to the 
relation. This fit is £(2695 A - V)/E(B - V) = 4.80 - 0.44 R(V). The nearly coincident dashed line is the 
fit in Figure [18] transformed by Equation (jnj). This fit is £'(2695 A - V)/E(B -V) = 4.73 - 0.42 R(V). 
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Fig. 20. — Plots of UV fit parameters (and derived quantities) vs. R(V). 1-a error bars are shown, based 
on our Monte Carlo simulations as in previous figures. The quantity A1250 is the difference between the 
observed value of E(X — V)/E(B — V) at 1250 A (8.0 /im" 1 ) and an extrapolation of the linear plus bump 
components of UV extinction, ft thus measures the strength of the FUV curvature and is computed by 
Af 250 = c 4 (8.0 — c 5 ) 2 (see Eq. [2]). Ab ump = tt 03/(27) is the area under the 2175 A bump. E\> ump = C3/7 2 
is the height of the 2175 A bump above the linear background extinction. 
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Fig. 21. — Plot of the slope of the UV linear component C2 vs. J?(F) _1 . 1-er error bars are shown, based on 
our Monte Carlo simulations as in previous figures. The solid line is a weighted linear fit which minimizes 
the residuals in the direction perpendicular to the fit. It is given by = 0.17 + 0.20 c 2 . 



